Surface spin flip probability of mesoscopic Ag wires 
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Spin relaxation in mesoscopic Ag wires in the diffusive transport regime is studied via nonlocal 
spin valve and Hanle effect measurements performed on permalloy/Ag lateral spin valves. The ratio 
between momentum and spin relaxation times is not constant at low temperatures. This can be 
explained with the Elliott- Yafet spin relaxation mechanism by considering the momentum surface 
relaxation time as being temperature dependent. We present a model to separately determine spin 
flip probabilities for phonon, impurity and surface scattering and find that the spin flip probability 
is highest for surface scattering. 
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Understanding how confinement influences physical 
properties is crucial for advancing nanotechnology 
Numerous studies have shown that when one or more 
dimensions of a structure become comparable to a char- 
acteristic length scale of a physical process in question 
(e.g., a mean free path for electron transport) even clas- 
sical boundary or surface effects can give rise to dramat- 
ically different behavior than that expected for the same 
bulk material. Examples include magnetoresistance in 
semiconductor nanostructures (negative vs. positive in 
the bulk) or thermal conductivities in Si nanowires 
(orders of magnitude reduction compared to bulk Si) Q . 
In contrast, confinement effects are less evident in metal- 
lic transport due to inherently short mean free paths but 
often manifest themselves in optical properties An 
important question to be addressed in spintronics Q 
is how does the size of a spin conductor or the surface 
conditions affect the transport of spin currents? Due to 
the relatively long spin diffusion length compared to the 
mean free path, confinement effects can be more pro- 
nounced in spin transport, even in metallic structures. 
So far, experiments performed with metallic lateral spin 
valve (LSV) structures @,0], where pure spin currents in 
a non-magnetic normal metal (N) are generated by dif- 
fusion of the non-equilibrium spin accumulation injected 
from a ferromagnet (F)Q, have focused mostly on deter- 
mining spin diffusion lengths l s and spin injection efficien- 
cies for various combinations of F/N materials, without 
quantifying contributions of different scattering mecha- 
nisms to the spin relaxation. In particular, to what ex- 
tent does confinement affect the spin relaxation time t s 
In this Letter we present a model, based on the 



Elliott- Yafet (EY) mechanism of spin relaxation [lfj, lll| , 
to separately quantify spin flip probabilities for phonon, 
impurity and surface scattering in mesoscopic metal wires 
in the diffusive transport regime. By studying spin trans- 
port in permalloy (Py)/Ag LSVs we find that the spin 
flip probability is highest for electron scattering from the 
Ag surface. Our model can also explain recent experi- 
mental results on temperature T [12j as well as thickness 
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FIG. 1: (Color online) (a) An SEM image of a Py/Ag LSV 
device adapted to show the nonlocal measurement configura- 
tion. Also shown are the directions of H\\ and H± applied in 
NLSV and Hanle effect measurements, respectively, (b) R n i 
vs. H\\ at 20 K. Corresponding M orientations of the Py elec- 
trodes are shown as blue arrows, while the total AR a signal 
is highlighted in red. (c) T dependencies of AR S and p. 



dependence of l a in mesoscopic Cu wires [13f . 

The Py/Ag LSV devices were fabricated on a SiN 
(100 nm)/Si substrate by e-beam lithography and shadow 
mask e-beam evaporation. A scanning electron micro- 
scope (SEM) image of a central region of the device is 
shown in Fig. 1(a). The two Py electrodes Pyl and Py2 
were both 25 nm thick and had widths of 130 and 80 nm 
respectively, while the bridging Ag wire was 260 nm wide 
and d = 80 nm thick. The center to center distance L 
between Py electrodes was 705 nm. Nonlocal spin valve 
(NLSV) and Hanle effect measurements were performed 
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by applying a dc current I = ± 0.3 mA from Pyl into 
the left part of the Ag wire and measuring the voltage 
V between Py2 and the right end of the Ag wire as a 
function of parallel H\\ and perpendicular H± magnetic 
fields, respectively [see Fig. 1(a)]. 

Figure 1 (b) shows the result of the NLS V measurement 
at 20 K. The dips in the nonlocal resistance, R n i = V/I, 
due to spin accumulation in Ag, upon switching the mag- 
netization M orientation of Pyl, are clearly observed. 
The magnitude of the difference between the R n i values 
measured for parallel and antiparallel M orientation of 
Py electrodes, AR S , is ~ 6.1 mO, which is a large sig- 
nal given that L — 705 nm [3, [H| . Also note that the 
values of R n i for parallel and antiparallel M orientation 
of the Py electrodes are almost perfectly symmetric with 
respect to zero, meaning that they are due to pure spin 
transport without parasitic ohmic signals [l6, 17 1. 

The large AR S values facilitate measurements of its 
T dependence by measuring R n i(T) for parallel and an- 
tiparallel remanent M orientations of the Py electrodes, 
respectively. Figure 1(c) shows the T dependence of 
AR„ = Rl] - 



I "t" 

Rtj and the corresponding conventional 



electrical resistivity p of the Ag wire. AR S is non- 
monotonic at low T despite the monotonic decrease of 
p [see Fig. 1(c)]. This behavior is consistently observed 
in all measured samples. In addition, similar behavior 
has also been observed in the case of Py/Cu LSVs 12 1 
and was attributed to the reduction of l s in Cu due to 
surface scattering. 

In order to determine whether a similar physical mech- 
anism also causes the behavior observed in our sam- 
ple, we turned to Hanle effect measurements, since they 
avoid variabilities between different samples, which are 
unavoidable in thickness-dependent studies. For Hanle 
measurements a combined effect of spin precession, relax- 
ation and dephasing leads to a characteristic dependence 
of AR S on H±, the Hanle resistance Rjj , given as: 18 1 



R$(H±) = tLEU J <p(t) CO s(w L t) exp dt, 

(1) 

for parallel M orientation of Py electrodes, and 
Rjj(H±) = —Rjj(H±) for antiparallel one. Here, 
ujl = gp,B^oH±/h is the Larmor frequency (g is 
the Lande factor, ps is the Bohr magneton, po is 
the magnetic permeability of free space) and V(t) — 
(l/V&nDt) exp(— L 2 /4Di) is the probability distribution 
of traveling times t of the injected spins from Pyl to Py2. 
Thus the Hanle effect measurements are used to deter- 
mine separately the injected spin polarization P and t s 
by fitting Rh vs. H± data to Eq. (1), if one knows the 
diffusion constant D, p, L, and the cross sectional area 
of the wire, A. 

Figure 2(a) shows representative Rh data, obtained 
at 40 K for Py electrodes prepared in parallel (red) and 
antiparallel (blue) M configurations. The difference be- 
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FIG. 2: (Color online) (a) Hanle signals measured at 40 K for 
parallel (i~t, red) and antiparallel (4-f, blue) M orientations 



of Py electrodes. Inset: ARh = Rjj ~ Rjj ■ Red line marks 
ARjj = for clarity, (b) sin 8 vs. poH± obtained from data 
shown in (a) by using Eq. (3). (c) Hanle signal (symbols) at 
40 K from the data shown in (a) by using Eq. (4). Best fit 
according to Eq. (1) shown as a red line. 



tween the two signals (not shown) exhibits an oscillating 
sign change as a function of H±, as expected from the 
spin precession. However, the behavior is more compli- 
cated than expected from Eq. (1). The striking feature of 
the data in Fig. 2(a) is the asymmetric shape of the two 
curves with respect to Rh — 0, which has not been ob- 
served in pre viously reported Hanle effect measurements 
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Mi- 

derstand this asymmetry, observed 
at all Ts, one has to take into account that in addition to 
precession and dephasing of the spin accumulation, the 
measured signal also depends on the orientation of the M 
of the Py electrodes with respect to the substrate plane. 
Namely, M inevitably tilts in the perpendicular direction 
due to application of H± . This decreases the fraction of 
precessing spin accumulation, and tends to restore the 
Rh signal to its initial value of AR s /2 for parallel M 
orientations. When this effect is taken into account, Rh 
can be expressed as: [l~8| 



R 



mm 

H 



(Ha 



= ±R f jj (H ± )cos 2 (e) + | (0)| sin 2 (<9), 

(2) 

with " + " and " — " signs corresponding to the an d 
4-t case, respectively. Here 9 is the angle between the 
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substrate plane and the direction of M. Based on Eq.(2): 
R%(H X ,6) + R%(H X ,6) 



sm z (9) 



2\R H (0)\ 



(3) 



Figure 2(b) shows the dependence of sin (9 on poH± 
obtained using Eq. (3) for the case shown in Fig. 2(a). 
For sufficiently low values of fioH± the dependence is lin- 
ear, but saturates above ~ ±0.5 T. This dependence is 
consistent with the Stoner-Wohlfarth model [2(J for co- 
herent M rotation with fields applied along a hard-axis 
direction. The slope of 2.7 T _1 for sin£> around zero 
field corresponds to a demagnetizing factor N = 0.37, 
taking fioM s = 1 T for Py. This value agrees reason- 
ably well with the literature one of N = 0.5 taking into 
account that the latter is defined for an infinitely long 
wire. Therefore, we conclude that the asymmetry in the 
measured Rjj curves arises from the tilting of the Py 
magnetizations. 

From the data, and Eq. (2), we can extract the Hanle 
signal: 



Rg(H x ) = \R„(0)\ 



Rh(H±,0) — R^ (H±,6) 



41'/ 



2\R h (0)\-Rh(H 



while 
Rg(H A 



(4) 

R l J(H ± ) = -R^(H±). Figure 2(c) shows 
) obtained using Eq. (4) for 40 K. The best fit ac- 
cording to Eq. (1) is shown as a red line. Only the region 
between ± 0.25 T is used since R^J(H±) is not well de- 
fined for higher fields due to the denominator in Eq. (4) 
being close to zero. The fit gives P — 0.207 and t s = 1 4.4 
ps, corresponding to l s = 564 nm. These values agree well 
with the ones reported in Ref. [3] (note that difference 
in l s scales with the difference in p) , but are considerably 
different than those reported in Ref. 15]. This suggests 
that determination of l s based on the transparent inter- 
face model used in Ref. [15( is not appropriate, due to 



the presence of an insulating oxide layer at the Py/Ag 
interface [2l[. We repeated the above procedure to de- 
termine P and t s in the range from 4.5 to 200 K, where 
the Hanle effect was observed. 

Figure 3(a) shows the T dependence of P. This depen- 
dence is monotonic and can be fitted to the Bloch model 
of thermally excited spin waves, i.e. P(T) = Po(l — nT 3 / 2 
as expected 
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In contrast, t s exhibits a maximum 
around ~40 K [see Fig. 3(b)] and then slightly decreases 
with decreasing T. This confirms that the observed non- 
monotonic T dependence of AR S [see Fig. 1(c)] is due to 
the reduction of r s , and hence l Sl at low T. Similar be- 
havior has been observed for Cu [12j . The dependence of 
spin relaxation rate l/r s on p is shown in Fig. 3(c). l/r s 
increases linearly with p above ~40 K, as expected from 
the EY mechanism of spin relaxation [l(| El, HH, but 
it exhibits a minimum around this T, and then slightly 
increases with decreasin g p. We point out that the sur- 
face spin relaxation 24j, |25j based on the Fuchs model 
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FIG. 3: (a) T dependence of P. The red line is fit to the 
Bloch model of thermally excited spin waves (see text), with 
P = 0.214 ± 0.003 and r\ = (4.0 ± 0.6) x 10~ 5 K~ 3/2 . (b) T 
dependence of t 3 . (c) Dependence of spin relaxation rate on 
resistivity of the Ag wire. 



26] with a T-independent surface momentum relaxation 
time cannot quantitatively explain the nonlinear de- 
pendence of 1/t s on p, let alone the upturn of l/r s at low 
T. The discrepancy can be resolved by invoking the con- 
cept of a T-dependent [27h2°| . Namely, when trans- 
port in the wire is diffusive, the , which is the average 
time it takes an electron to diffuse over the distance d is 
given as r/ = 7c? 2 /Db- Here, 7 is the averaging coeffi- 
cient and Db = (1/2>)v 2 f t® is the T-dependent bulk dif- 
fusion constant, with vf = 1-39 x 10 6 m/s and being 
the electron Fermi velocity and the bulk momentum re- 
laxation time, respectively, . The latter is determined by 
scattering within the bulk of the wire and can be defined 
as (jf)- 1 = (7f h ) -1 + (T l e mp )-\ where rf and r e imp are 
momentum relaxation times for electron scattering from 
phonons and impurities (including the grain boundaries), 
respectively. Thus, the total momentum relaxation time 
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Associating to each scattering process its corresponding 
spin flip probability, i.e., e p h, timp and es, and follow- 
ing the EY proportionality between momentum and spin 
relaxation times 1231 . we find 
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(6) 
(7) 



Equation (7) naturally explains the nonlinearity in l/r s 
vs. p (since within this model depends nonlinearly on 
p) and, even more significantly, it describes the upturn 
in l/r s at low T, since the first and second terms in 
Eq. (7) have different dependencies on p. Also, based on 
Eq. (7), one can determine e p h, £imp and es by fitting 
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FIG. 4: (Color online) Dependence of spin relaxation rate in 
the Ag wire on the bulk momentum relaxation rate (symbols) . 
The fit to the data according to Eq. (6) is shown as a red line. 



the 1/t s vs. dependence, which can be obtained by 

performing T-dependent spin transport measurements, 
such as is shown in Fig. 3. 

Figure 4 shows the plot of l/r s vs. (black cir- 

cles), where is obtained using Eq. (5) as = (1 — 
a/I -4aT e 2 )/(2aT e ), with a = {v F /d) 2 /3j = 179.3 ps~ 2 
was obtained by determining 7 = 0.563 using Dingle's ef- 
fective mean free path model for completely diffuse elec- 
tron surface scattering (jjjj, and r e = m/(ne 2 p) was de- 
termined from the measured p using n = 5.85 x 10 28 m~ 3 . 
The best fit according to Eq. (7) is shown as a red line. 
We find e ph = (7.5 ± 1.3) x 1CT 3 , e s = (1.7 ± 0.4) x 10~ 2 
and Cimp = (—5.0 ± 5.2) x 10~ 3 . These values show 
that the spin flip probability is highest for surface scat- 
tering and weakest for scattering from impurities (< 
2 x 10 _4 )[32J. Also note that the value for e p h agrees 
well with a previously reported value of e p h = 2 x 10~ 3 
for Cu wires (~ 4x smaller value, consistent with Cu 
being a lighter element) 0, [HI as well as the value of 
e p h = 2.9 x 10 -3 obtained from electron spin resonance 
on bulk Ag [ilj]. 

Based on this model, we can also predict the thickness 
dependence of l s . Multiplying Eq. (7) by 1/D one finds 



I, = 



^Jad 2 + 



(8) 



r tmp 



where in the low T limit, when rf 
e mp /{DTl m P) and f3 = (v F e S T™P)/(3jD). In the high 
T limit, where phonon scattering is non-negligible, a = 
eimp/iDr^P) + e ph /(DTP h ), and /? is the same. These 
relations in principle can be used to determine e p h, timp 
and es by fitting the thickness dependence of the spin 
diffusion length obtained for low and high T. However, 
such measurements require several samples, which intro- 
duce additional experimental uncertainties. 

In conclusion, we have studied the spin relaxation in 
a mesoscopic Ag wire in the diffusive transport regime, 



and observed a nonlinear dependence of the spin relax- 
ation rate on resistivity. This observation cannot be ex- 
plained quantitatively with the Elliott- Yafet mechanism 
of spin relaxation by a conventional approach, which con- 
siders surface spin relaxation as being temperature inde- 
pendent. We present a model that explains these ob- 
servations with the Elliott- Yafet mechanism by adding a 
temperature dependence to the surface relaxation. This 
enables us to quantify spin flip probabilities for phonon, 
impurity and surface scattering respectively. We find 
that the spin flip probability of the Ag wire is strongest 
for electron scattering from surfaces. 
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